First-principles study of ground state properties of ZrH2 
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Abstract 

Structural, mechanical, electronic, and thermodynamic properties of fluorite and tetragonal 
| ; phases of Z r H 2 „ syst e maMy , udl e d by e mpl o ymg the den* f— . tW y within g en- 

eralized gradient approximation. The existence of the bistable structure for ZrH2 is mainly due to 
the tetragonal distortions. And our calculated lattice constants for the stable face-centered tetrag- 
onal (fct) phase with c/a=0.885 are consistent well with experiments. Through calculating elastic 

constants, the mechanically unstable characters of face-centered cubic (fee) phase and fct structure 

> '. 

with c/a=l.lll are predicted. As for fctO.885 structure, our calculated elastic constants explicitly 

CO 

in 



Ctf 



indicate that it is mechanically stable. Elastic moduli, Poisson's ratio, and Debye temperature are 



derived from elastic constants. After analyzing total and partial densities of states and valence 

o : 

electron charge distribution, we conclude that the Zr— H bonds in ZrH2 exhibit weak covalent fea- 
ture. But the ionic property is evident with about 1.5 electrons transferring from each Zr atom to 

•i-H , 

H. Phonon spectrum results indicate that fctO.885 and fctl.lll structures are dynamically stable, 

5-H 



while the fee structure is unstable. 

PACS numbers: 61.50.Ah, 62.20.Dc, 71.15.Mb, 63.20.dk 
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I. INTRODUCTION 

Many transition metals react readily with hydrogen to form stable metal hydrides |l|. 
The common metal hydrides are technologically attractive materials due to their ability 
to store high densities of hydrogen safely. In the case of zirconium hydrides, another ap- 
pealing interest arises in nuclear industry, essentially used as a neutron moderator [2J and 
fuel rod cladding materials |3| in nuclear reactors. Additionally, in the fusion technology 
applications, zirconium hydrides are unavoidably precipitated in the excess of the hydrogen 
solubility, which may adversely affect several mechanical and thermal properties in some 
crystal orientations [4] involving different stable and met ast able hydride phases which may 
lead to a significant embrittlement. Furthermore, zirconium hydrides are also served as 
a part of the promising new type of actinide hydride fuels, such as U-(Th-Np-Am)-Zr-H 
[5j. Consequently, great attentions are needed to focus on the basic scientific research of 
zirconium hydrides. 

Since 1952, there has occurred in literature large amounts of experimental investigations 
on zirconium hydrides. The electronic structures of ZrH x . (1.63<x<1.94) was studied using 
photoelectron spectroscopy and synchrotron radiation [6J, where the Jahn- Teller effect is 
demonstrated by the changes of electrons occupation at the Fermi level. An assessed H-Zr 

n. 

phase diagram was successfully constructed by Zuzek et al. [7[ in 1990. And the transition 
between the cubic and the tetragonal phases has been suggested by Ducastelle et al. [8fl to be 
of Jahn- Teller type too. From theoretical point of view, the phase transition between cubic 
and tetragonal phases has been investigated by Switendick J9(] , where the non-self-consistent 
augmented plane wave (APW) method was employed. Most importantly, the study of the 
electronic structure and energetics for the tetragonal distortion in ZrH2 has been extensively 
performed. Using the pseudopotential method with both local density approximation (LDA) 
and general gradient approximation (GGA), Ackland [10] obtained a double minimum in the 
energy surface, one minimum for c/a<l and the other for c/a>l with an energy barrier in 
between of roughly 0.1 eV. But the ground state for c/a>l was in clear disagreement with 



the experimental data llMlq]. Recently, Wolf and Herzig [17] and Quijano et al. [18] studied 
the accurate total energy and the electronic structure characteristic of ZrH_2 by using LDA 
and GGA, respectively, within the full potential linear augmented plane wave (FP-LAPW) 
method. Both of these works obtained the correct ground state for c/a<l and made a clear 



explanation of the cubic to tetragonal distortion. 

However, despite the abundant theoretical and experimental research on zirconium hy- 
drides, relatively little is known regarding their chemical bonding nature, mechanical prop- 
erties, and the phonon dispersions. Until now, the elastic properties, which relate to various 
fundamental solid state properties such as interatomic potentials, equation of states, phonon 
spectra, and thermodynamical properties are nearly unkown for ZrH^. Moreover, although 
the electronic properties as well as the chemical bonding in zirconium dihydride have been 
calculated 17], the study of the bonding nature of Zr-H involving its mixed ionic/covalent 
character is still lacking. As a consequence, these facts inhibit deep understanding of the 
zirconium dihydride. Motivated by these observations, in this paper, we present a first- 
principles study by calculating the structural, electronic, mechanical, and thermodynamical 
properties of zirconium dihydride. Our calculated results show that the bistable structure 
of ZrH 2 is the fct phase at c/a=0.885 and c/a=l.lll, respectively. The mechanical and dy- 
namical stability of fct and fee phases of ZrH2 are carefully analyzed. In addition, through 
Bader analysis |l9l . [20J we find that about 1.5 electrons transfer from each Zr atom to H 



atom for ZrH.2- 



II. COMPUTATIONAL METHOD 



Our total energy calculations are self-consistently carried out using density functional 



theory (DFT) as implemented in Vienna ab initio simulation package (VASP) 



2l|, which is 



based on pseudopotentials and plane wave basis functions. All electron projected augmented 



wave (PAW) method of Blochl |22j is applied in VASP with the frozen core approximation. 
The generalized gradient approximation (GGA) introduced by Perdew, Burke, and Ernz- 
erhof (PBE) [23] is employed to evaluate the electron exchange and correlation potential. 
Zirconium 4s 2 4p 6 4d 2 5s 2 and hydrogen Is 1 electrons are treated as valence electrons. The 
integration over the Brillouin Zone (BZ) is performed with a grid of special k point-mesh de- 



termined according to the Monkhorst-Pack scheme [24|. After convergence test, 12x12x12 k 
point-mesh is chosen to make sure the total energy difference less than 1 meV per primitive 
cell. When exploring the electronic structure, a finer k point-mesh 24x24x24 is preferred 
to generate a high quality charge density. The expansion of the valence electron plane wave 
functions can be truncated at the cutoff energy of 650 eV, which has been tested to be fully 




FIG. 1: Tetragonal unit cell in space group IA/mmm (a) and cubic unit cell for ZrH2 in space 
group Fm3m (b) with larger cyan spheres for Zr atoms and the smaller white H. 



converged with respect to total energy. Relaxation 



out according to the Methfessel-Paxton scheme 25J with a smearing width of 0.1 eV. The 



procedures at ground state are carried 



geometry relaxation is considered to be completed when the Hellmann-Feynman forces on 
all atoms were less than 0.001 eV/A. And the accurate total energy calculations are per- 
formed by means of the linear tetrahedron method with the Blochl's corrections 26| . The 
self-consistent convergence of the total energy calculation is set to 10 -6 eV. 



III. RESULTS 



A. Atomic structure and mechanical properties 



At ambient condition, the stable zirconium dihydride crystallizes in a fct structure with 
space group I4/mmm (No. 139). In its unit cell, there are two ZrH2 formula units with Zr 



and H atoms in 2a:(0,0,0) and 4d:(0,|,j) sites, respectively [see Fig. 1(a)]. Each Zr atom 
is surrounded by eight H atoms forming a tetragonal and each H connects with four Zr 
atoms to build a tetrahedron. The present optimized equilibrium lattice parameters (a and 
c) obtained by fitting the energy-volume data in the third-order Birch- Murnaghan equation 
of states (EOS) 27| are 5.030 A and 4.414 A (see Table I), in good agreement with the 
experimental values. In addition, the fee fluorite type structure with space group Fm3m 
(No. 225) is considered as the metastable phase for ZrH2. The cubic unit cell is composed 
of four ZrH 2 formula units with the Zr and H atoms in 4a:(0,0,0) and 8c:(|,|,|) sites, 
respectively [see Fig. 1(b)]. Here, each Zr atom is surrounded by eight H atoms forming a 
cube and each H atom connects with four Zr atoms to build a tetrahedron. Our optimized 
equilibrium lattice constant a for fee ZrH2 is 4.823 A. This value, although no experimental 



TABLE I: Calculated lattice constants, bulk modulus Bq, and pressure derivative of the bulk 
modulus Bq at ground state from EOS fitting. For comparison, available experimental values and 
other theoretical results are also listed. 



Compounds Property Present work Previous calculation 



Experiment 



fee ZrH2 



fct ZrH 2 



ao (A) 


4.823 


B (GPa) 


133 


B'o 


4.01 


a (A) 


5.030 


co (A) 


4.414 


B (GPa) 


127 


B' 


3.62 



4.817 a , 4.804 & 
136 a ,152 fe 

5.021 a , 5.008 fe , 5.000 c 4.985**, 4.975 e , 4.982^ 
4.432 a , 4.419 & , 4.450 c 4.430 d , 4.447 e , 4.449 / 



Reference [l8j] , b Reference 17|] , c Reference 



m. 
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FIG. 2: Total energy as a function of the c/a ratio in the fct phase for ZrH\2. 

structural values are available for comparison, is in good agreement with previous theoretical 
work in the literature (see Table I). 

In present study, the stability of ZrH2 has been investigated by calculating the total 
energy as a function of c/a at the optimized constant volume of fee structure, as illustrated 
in Fig. 2. The energy curve displays two local minima, one with c/a<l (0.885) and the other 
with c/a>l (1.111). The minimum at c/a=0.885 has the lowest energy and corresponds to 
fct phase observed experimentally. The fet-fee structure energy barrier obtained from this 
calculation is 0.1048 eV, in good agreement with the Ackland's results 10] but much larger 



than the FP-LAPW calculation results 
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TABLE II: The strain combinations in the strain tensor [Eq. (1)] to calculate the elastic constants 
of cubic and tetragonal ZrH 2 . 



Phase Strain Parameters (unlisted ej=0) 



AE/V in 0(5 2 ) 



fee 



fct 



ei=e 2 =e 3 =5 




(|C U + 3C 12 ),5 2 


ei=e2=S 




(Cii + C 12 )5 2 


e 4 =e 5 =e 6 =5 




2C440 


ei = e 2 =e 3 =<5 


(c 


li + C 12 + 2C 13 + ±C 33 )<5 2 


ei=e 3 =5 




±(Cn+2C 13 + C 33 )<5 2 


ei=e 2 =c) 




(Cn + C 12 )5 2 


e 3 =5 




^ 33 <5 2 


e 4 =e 5 =5 




C M 5 2 


e 6 =5 




^CqqS 



Elastic constants can measure the resistance and mechanical features of crystal to external 
stress or pressure, thus may evaluate the stability of crystals against elastic deformation. 
For small strain e, Hooke's law is valid. Therefore, we can enforce it onto the equilibrium 
lattice, determine the resulting change in the total energy, and from this information deduce 
;he elastic constants. The crystal total energy E(V,e) can be expanded as a Taylor series 

6 



E(V, e) = E(V , o) + v J2 n* + y E C ^ e i + °({ e i})> W 

where E(Vq, 0) is the energy of the unstrained system with the equilibrium volume Vo, e is 
the strain tensor which has matrix elements Sij (i,j—l, 2, and 3) defined by 



I P , £6 £5 \ 
C1 2 2 

9 °2 9 



\f f C3 



(2) 



and Cy are the elastic constants. Such a strain transforms the three lattice vectors defining 
the unstrained Bravais lattice {a^, k—1, 2, and3} to the strained vectors {a' fc } [29| as defined 
by 



fl + e)a*, 



(3) 



TABLE III: Calculated elastic constants in units of GPa for ZrH 2 at zero pressure. 



Phase C\\ Ci 2 C44 C13 C33 Cqq 

fee 82.6 159.7 -19.5 

fctO.885 165.6 140.9 30.5 106.8 145.5 60.6 

fctl.lll 125.7 145.5 30.9 115.0 190.6 42.0 

where I is a unit 3x3 matrix. Each lattice vector a& or a' fc is a 3 x 1 matrix. 

For the cubic phase, there are three independent elastic constants, i.e., Cn, C i2 , and C 44 , 
which can be calculated through a proper choice of the set of strains {ej,z = 1, ...,6} (see 
Table II). As for the fct phase, the six independent elastic constants, i.e., Cn, Ci 2 , C13, 
C33, C44, and 6*66, can be obtained from six different strains listed in Table II. To avoid the 
influence of high order terms on the estimated elastic constants, we have used very small 
strains, i.e. within ±1.5%. Thus, in our calculations the strain amplitude 5 is varied in 
steps of 0.006 from 5=— 0.036 to 0.036. And the total energies E(V, S) at these strain steps 
are calculated and fitted through the strains with the corresponding parabolic equations of 
AE/Vo to yield the required second-order elastic constants. 

Our calculated elastic constants for fee ZrH 2 and fct ZrH2 at c/a=0.885 and c/a— 1.111 
are collected in Table III. Obviously, fct ZrH 2 at c/a=0.885 is mechanically stable due to 



the fact that its elastic constants satisfy the following mechanical stability criteria 30( of 
tetragonal structure: 

Cn > 0, C 33 > 0, C 44 > 0, C 66 > 0, 
(C n - C 12 ) > 0, (C n + C 33 - 2C 13 ) > 0, 

[2(C 11 + C 12 ) + C 33 + 4C 13 ]>0. (4) 

Nevertheless, the fct phase at c/a=l.lll is mechanically unstable in light of the fact that 
Cn is smaller than Ci 2 . As for the fee structure of ZrH 2 , it is also mechanically unstable 
because of the fact that its elastic constants cannot meet the following mechanical stability 



criteria [30f| of cubic structure: 

C n > 0,C 44 > 0,C n > |C 12 |, (C n + 2C 12 ) > 0. (5) 

After obtaining elastic constants, bulk modulus B and shear modulus G for fct phase 
of ZrH 2 at c/a=0.885 can be calculated from the Voigt-Reuss-Hill (VRH) approximations 



TABLE IV: Calculated bulk modulus, shear modulus, Young's modulus, Poisson's ratio, density, 
transverse sound velocity, longitudinal sound velocity, average sound velocity, and Debye temper- 
ature of the fct phase of ZrH2 at c/a=0.885. 



5(GPa) G(GPa) 


£(GPa) 


V 


p(g/cm 3 ) 


vt{km/s) 


vi (km/s) 


« m (km/s) D (K) 


130 29 


80 


0.397 


5.5201 


2.2823 


5.5197 


2.5825 364.9 



3lN33|. Then the Young's modulus E and Poisson's ratio v are calculated through E = 
9BG/(3B + G) and v = (3B - 2G)/[2(3B + G)]. The Debye temperature (9 D ) is obtained 
using the relation 34] 

1/3 



e 



D 



h 
k B 



3n 

Ait 



m 



(6) 



where h and fc# are Planck's and Boltzmann's constants, respectively, Na is Avogadro's 
number, p is the density, M is the molecular weight, n is the number of atoms in the 
molecule, and v m is the average wave velocity. The average wave velocity in the polycrys- 
talline materials is approximately given as 

1 1 -1/3 



-( — —' 

3V + ^f j 



(7) 



where v t =\/G j ' p and vi=y/(3B + AG)/3p are the transverse and longitudinal elastic wave 
velocity of the polycrystalline materials, respectively. 

All the calculated results for fctO.885 are collected in Table IV. Results of fee and fct 1.1 11 
can not be obtained due to their mechanically unstable nature. Note that we have also 
calculated the bulk modulus B by EOS fitting. The derived bulk modulus for fctO.885 turns 
out to be very close to the one obtained from the above VRH approximation, which again 
indicates that our calculations are consistant and reliable. It is well known that the shear 
modulus G represents the resistance to plastic deformation, while the bulk modulus B can 
represent the resistance to fracture. A high (low) B/G value is responsible for the ductility 
(brittleness) of polycrystalline materials. The critical value to separate ductile and brittle 
materials is about 1.75. Using the calculated values of bulk modulus B and shear modulus 
G for fctO.885, the B/G value of 4.48 can be obtained. For element Zr, the B/G value of 
2.63 can be derived from the elastic data in Ref. 35J. Therefore, fct phase of ZrH 2 is rather 



ductile and its ductility is more predominant than that of element Zr. As for the Poisson's 
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FIG. 3: Total and orbital-resolved local densities of states for (a) fee, (b) fctO.885, and (c) fctl.lll 
structures of ZrH2. The Fermi energy level is set at zero. 

ratio, the value is well within the range from 0.25 to 0.45 for typical metals. Unfortunately, 
no reliable experimental and theoretical results concerning on the mechanical properties in 
the literature are available for comparison. We hope that our calculated elastic constants 
and elastic moduli can be illustrative in the realistic application of the mechanical data for 
ZrH 2 . 



B. Electronic structure and charge distribution 

Basically, all the macroscopical properties of materials, such as hardness, elasticity, and 
conductivity, originate from their electronic structure properties as well as the nature of the 
chemical bonding. Therefore, it is necessary to perform the electronic structure analysis of 
ZrH^. The calculated total densities of states (DOS) and the orbital-resolved partial densities 
of states (PDOS) of fee ZrH_2, fct Zrri2 at c/a=0.885 and c/a— 1.111, are shown in Fig. 3. 
Wholly, the occupation properties of ZrH 2 in fct and fee phases are similar. The pseudogap 
below the Fermi level, located at around —3.5 eV, for all three systems can be observed. The 
energy region below the pseudogap has predominant feature of hybridization for H Is orbital 
and Zr 4d and Zr 4p orbitals. States in this region have critical contribution to the Zr-H and 
H-H bonding state. However, from the pseudogap to the Fermi level, His states contribute 
very little to the total DOS, while the Zr 4d states dominate in this energy range. The fact 
that the DOS occupation has obvious sharp peak at the Fermi level for fee ZrH 2 implies that 
the cubic phase should be unstable. In contrast, a strong reduction in the DOS at the Fermi 
level is observed in the two fct structures. This evident difference between fee and fct phase 




FIG. 4: Valence charge density of (a) fee phase in (110) plane, (b) fctO.885 in (100) plane, and (c) 
fctl.lll in (100) plane. The contour lines are drawn from 0.0 to 1.0 at 0.05 e/A 3 intervals. 

in the vicinity of the Fermi level is caused by the splitting of the band in this region. As 
a consequence, the Fermi level in the two tetragonal phases is not located at a peak of the 
DOS but at a local minimum, which leads to a reduction in the electronic contribution to the 
total energy. This behavior has been referred to as a Jahn- Teller mechanism, as indicated 
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18| 



by previous experiments and theoretical works p, ll2|, 

To analyze the ionic/covalent character of zirconium dihydride, in the following we will 
carefully investigate the valence charge density distribution. The calculated valence charge 
density maps of the fct ZrH2 at c/a=0.885 and c/a=l.lll in (100) plane and fee Zrli~2 in 
(110) plane are plotted in Fig. 4. Obviously, the charge densities around Zr and H ions are 
all near spherical distribution with slightly deformed towards the direction to their nearest 
neighboring atoms. Covalent bridges, which represent the bonding nature of Zr— H bonds, 
is clearly shown in Fig. 4. In Fig. 5, we plot the line charge density distribution along 
the nearest Zr— H bonds. Clearly, three line charge density curves vary in almost the same 
way along the Zr— H bonds. One can find a minimum value of charge density for each bond 
at around the bridge locus (indicated by the arrow in Fig. 5). These minimum values are 
listed in Table V. Although these values are much smaller than 0.7 e/A for Si covalent bond, 
they are prominently higher than 0.05 e/A for Na— CI bond in typical ionic crystal NaCl. 
Therefore, there are weak but clear weights in Zr— H bonds in both fee and fct phases of 
ZrH 2 . 

In addition, we have performed the Bader analysis |19l . |20j] for the three typical ZrH2 
cells. The charge (Qb) enclosed within the Bader volume (Vb) is a good approximation 
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FIG. 5: The line charge density distribution between Zr atom and the nearest neighbor H atom 
for fee, fctO.885, and fctl. Ill structures of ZrH2. 

to the total electronic charge of an atom. In this work, the default charge density grids 
for one primitive cell are 56x56x56, 56x56x56, 60x60x60 for fee, fctO.885, and fctl. Ill 
cells, respectively. To check the precision, the charge density distributions are calculated 
with a series of n times finer grids (n=2,3,4,5,6,7). The deviation of the effective charge 
between the six and the seven times finer grids is less than 0.2%. Thus we perform the 
charge density calculations using the seven times finer grid (392x392x392, 392x392x392, 
420x420x420) for fee, fctO.885, and fctl. Ill cells, respectively. The calculated results are 
presented in Table V. Note that although we have included the core charge in charge density 
calculations, only the valence charge are listed since we do not expect variations as far as 
the trends are concerned. From Table V, we find that the Bader charges and volumes for 
fee and fct ZrH2 are almost equal to each other. This shows similar ionic character, through 
a flux of charge (about 1.5 electrons for each Zr atom) from cations towards anions, for 
zirconium dihydride in their stable phase and metastable phases. The slight difference in 
charge distribution is mainly due to the structure distortion. Thus the ionicity of Zr-H 
bonds for ZrH2 is also evident either for its fee structure or fct structures. 

C. Phonon dispersion curves 

Phonon frequencies of the crystalline structure is one of the basic aspects when considering 
the phase stability, phase transformations, and thermodynamics of crystalline materials. 
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TABLE V: Calculated charges Qb and volumes Vs(A 3 ) according to Bader partitioning as well as 
the Zr— H distance (A) and relevant minimum values of charge densities (e/A 3 ) along the Zr— H 
bonds for ZrH2. 

Compounds QbCZit) <5b(H) Fs(Zr) Vb(H) Zr-H distance Charge density m j n 

fee 10.484 1.758 14.809 6.618 

fctO.885 10.512 1.744 14.953 6.544 

fctl.lll 10.514 1.743 14.954 6.545 



2.088 


0.320 


2.095 


0.316 


2.094 


0.317 



Employing the Hellmann-Feynman theorem J36|, |37J and the direct method 38(, we have 
calculated the phonon curves along some high symmetry directions in the BZ, together with 
the phonon DOS. For the phonon dispersion calculation, we use the 2x2x2 fee (fct) supercell 
containing 96 (48) atoms for fee (fct) ZrH 2 and the 4x4x4 Monkhorst-Pack /c-point mesh 
for the BZ integration. In order to calculate the Hellmann-Feynman forces, we displace four 
and eight atoms, respectively, for fee and fct ZrH 2 from their equilibrium positions and the 
amplitude of all the displacements is 0.03A. The calculated phonon dispersion curves along 
the T-X-K-T-L-X-W-L directions for fee ZrH 2 and along the T-N-P-X-T-Z directions for 
fct ZrH 2 are displayed in Fig. 6. For both fee and fct ZrH 2 , there are only three atoms in 
their primitive cell. Therefore, nine phonon modes exist in the dispersion relations. Due to 
the fact that zirconium is much heavier than hydrogen, the vibration frequency of zirconium 
atoms is apparently lower than that of hydrogen atoms. As a result, evident gap between the 
optic modes and the acoustic branches exits and the phonon DOS of fee (fct) ZrH 2 can be 
viewed as two parts. One is the part lower than 7.4 (7.1) THz, where the main contribution 
comes from the zirconium sublattice, while the other part range from 30.7 to 37.9 (29.6 to 
38.5) THz is dominated by the dynamics of the light hydrogen atoms. As shown by Fig. 
6(b) and 6(c), all frequencies are positive for the two fct structures. This assures that the 
two fct phases are both dynamically stable. In contrast, one can see from Fig. 6(a) that 
the fee phase of ZrH 2 is dynamically unstable. The transverse acoustic (TA) mode close 
to T point becomes imaginary along the T-K and T-L directions. The minimum of the TA 
branch locates at the T-K direction. Therefore, the fee to fct phase transition probably can 
occur along the [101] direction. This kind of dynamical instability of fee phase compared to 
the stable fct phase is well consistent with our above-discussed mechanical stability analysis 
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FIG. 6: Calculated phonon dispersion curves along the high-symmetry directions (left panel) and 
the corresponding PDOS (right panel) at ground state for (a) fee ZrH2, (b) fct ZrH2 at c/a=0.885, 
and (c) fct ZrH2 at c/o=l.lll. 

of ZrH 2 . 

IV. CONCLUSION 

In summary, we have investigated the structural, mechanical, electronic, and thermody- 
namic properties of ZtB.2 in its stable and metastable phases by means of first-principles 
DFT-GGA method. Our optimized structural parameters are well consistent with previous 
experiments and theoretical calculations. Tetragonal distortion leads to bistable structures 
with c/a>l and c/a<l, the latter of which corresponds to the experimentally observed 
low-temperature phase. Our total energy calculations illustrate that the most stable phase 
of ZrH 2 is the fct structure with c/a=0.885. Elastic constants, various moduli, Poisson's 
ratio, and Debye temperature are calculated for fee and two fct phases of ZrldV Mechani- 
cally unstable nature for fee and fctl.lll structures are predicted. For fctO.885, mechanical 
stability and dynamical stability are firstly predicted by our elastic constants results and 
phonon dispersion calculation. The occupation characters of electronic orbital also accord 
well with experiments. Through Bader analysis, we have found that the Zr— H bonds ex- 
hibit weak covalent character while the ionic property is dominant with about 1.5 electrons 
transferring from each Zr atom to H. Our calculated phonon curves of fee ZrH2 have shown 
that the TA mode becomes imaginary close to T point along the [101] direction. 
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